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During the past five years there has teen a significant increase of 
interest in the problem of radiant heat transfer to absorbing gases. 
Recent activity in this area has been motivated by the advent of heat 
transfer problems arising in space vehicle reentry, MHD energy conver- 
sion, and energy transport in a gaseous nuclear reactor. 
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Hottel was one of the earliest workers in this field, having inves- 
tigated radiant heat transfer from furnace gases as early as 1927. 
Although other workers have recently become active in the field by ob- 
taining analytical solutions for heat transfer and temperature distribu- 
tions ^1), (2j, (5), Hot lei r s analysis ( 4) of the heat transfer problem 
remains probably the most realistic since it is applicable to real 
(nongray) gases and may be applied to geometries of almost arbitrary 
shape. 

There has been very little done so far in obtaining generalized 
results for radiation to a flowing gas or for combined radiation and 
conduction. Adrianov, (3), investigated radiation heat transfer to a 
flowing gas, but some of the assumptions in his analysis render his 
results valid only for the case of weak absorbtion in the gas. Viskanta, 
has presented a rather complete analytical treatment of combined 
radiation and conduction between two infinite parallel plates separated,,* 
by an absorbing gas, but only a few specific results were presented. 

It is the purpose of this investigation to expand on the work done 
by both Adrianov and Viskanta which is cited above. The method used in 
the present analysis is a modification of Hottel T s zoning technique ( 4 ). 
Results are obtained for a gray gas in a rectangular channel formed by 
two black, parallel flat plates of finite length, and are presented in 
generalized form for the case of radiation to flowing gases and also for 
combined radiation and conduction. 


ANALYSIS 

A two-dimensional analysis of radiation heat transfer is made for 
a gray gas enclosed in a rectangular channel formed by two black parallel 
flat plates of finite length and infinite width. A sketch of this con- 
figuration is shown in figure 1. A rigorous treatment of this problem 
requires the solution of the following two dimension integro-differential 
equation which represents the heat balance on an infinitesimal volume, 
dV, at any point in the channel. 
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absorbtion coefficient of the gas 
energy emitted per unit volume at r = r Q 

rate of enthalpy increase of the flowing gas at r - r 
net conduction heat transfer per unit volume 


integrals on radiation absorbed from emission due to surrounding gas 
the right volume and radiation absorbed from emission of flat plate 

hand side surfaces , respectively 

If the conduction and convection terms are removed in equation (l), its 
solution is greatly simplified since the equation will then be linear in 
the emissive power oT 4 . If in addition the length of the channel also 

b e co mes - Inf mite j - the - integral s canr be par tially - evaluated in el osed fo r m 

in terms of exponential -integral functions, as was done by Usiskin (1). 

In the absence of these simplifying restrictions, the only feasible 
way to solve this problem appears to be application of a method similar to 
Hottel's where the two-dimensional integral-differential equation (l) is 
approximated by a system of algebraic equations. This is done by first 
dividing the region between the plates up into a 10x10 array of rectangular 
zones, 100 zones in all. The gas temperature variation within each of 
these zones may be approximated by a two-dimensional linear function. 

Heat balance equations in the form of equation (l) are now written for 
infinitesimal volume elements located at the center of each of 100 zones. 
The integral on the right of equation (l) is approximated by a finite 
summation over all the zones and surfaces. The derivatives are replaced 
by numerical approximations obtained from fitting a paraboloid through 
the temperatures of the zone in question and those of its adjacent 
neighbors in both directions. Each of the 100 zones is labeled by a dual 
subscript (i,j), i is the position along the length of the channel and j 
represents the position in the transverse direction. The resulting system 
of nonlinear algebraic equations, which is given below, are then solved 
for the 100 values of Tj_j, the gas temperatures at the center of the 
zones, for a given set of boundary conditions (plate and end temperatures). 
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The conditions of the two ends of the channel are represented by assuming 
that the ends are porous black plugs having the same temperature as the 
effective gas temperature at the ends. The equations to be solved are thus: 
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oTg(m,n)g(r(m,n) - r(i,j)) i=l,10; j=l,10 (2) 


The principal difference between the present approach and that of Hottel's 
is that here the heat balance is taken on an infinitesimal volume of the 
center of the zone, rather than on the entire zone. This approach results 
in a large reduction in the effort required to compute the "exchange 
factors" f(r) and g(r) . Since the system of 100 equations (2) is gen- 
erally nonlinear due to the presence of convection and conduction terms, 
it was solved using the Newton -Raphson iterative method for nonlinear 
algebraic equations. This method rapidly converges and the correct solu- 
tion is obtained after only a few iterations. 

RESULTS AND DISCUSSION 

The system of equations described by equation (2) was^ solved _on__ an 

IBM 7090 computer for a gray gas of constant absorptivity enclosed in a 
black walled rectangular channel of aspect ratio, L/D, of ten and for a 
range of optical thickness t , (0.1 ^ t q ^ 6), t q = kD, where k is the 
gas absorptivity and D is the plate spacing. Cases were run for either 
combined convection and radiation or combined conduction and radiation, for 
a wide range of gas flow rates and thermal conductivities. The walls of the 
channel were considered to be at constant temperatures. To establish a 
reference several cases were run for radiation only, without flow or 
conduction. 

Figure 2 illustrates the heat transfer between the two plates for the 
case of a stagnant nonconducting gas. Since in this case equation (l) is 
linear in the emissive power, oT 4 , the heat transfer results may be normal- 
ized by dividing by the difference in the emissive powers of the two plates. 
The results given in figure 2 compare favorably with those obtained by 
Usiskin, (l), for infinite parallel flat plates. The agreement is not 
perfect because of end effects due to the finite length of the channel. 

In order to discuss the effect of conduction on radiation in a stagnant 
gas it is desirable to present the results in some type of nondimensional 
parametric form to obtain complete generality. Equation (l) may be 
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nondimen s i onali z ed by dividing through by kcrT^ where T* is the tempera- 
ture of the hotter plate. It can also be shown that the integrals on the 
right, for a given temperature field, are a function only of t q , and the 
channel aspect ratio L/D. On the basis of the above observations it can be 
shown that the solution of equation (1) is a unique function of the boundary 
conditions (ratio of cold plate and end temperatures to T*), and of the 


following set of dimensionless parameters 

A/D 


ft 

vcrT. 


T o> Tj)- 


For the sake of 


brevity the parameter will be designated N^. Figure 3 shows the 


effect of combined conduction and radiation on the transverse temperature 
profile for the cases of both weak gas radiation (t = 0.2), and a relatively 
opaque gas (t = 3. C) . In both cases, the temperature gradient and therefore 


the heat conduction at the cold wall were substantially higher than those for 
pure conduction (linear temperature profile). At the hot wall the temperature 
gradient may be either greater or less than that for pure conduction, depend- 
ing upon the value of and the opacity of the gas. 


Another effect of combined radiation and conduction is that the total 
heat transferred is greater than the sum of the heat transfer for radiation 
and conduction, taken separately. This effect is shown in figure 4 where 
the ratio of combined heat transfer to the sum of both components taken 
separately is plotted versus r Q for various values of N^. This augmenta- 
tion is due to interaction between radiation and conduction resulting from 
the nonlinear nature of the combined process. There is no augmentation for 
the case t q = 0 since the gas is then not involved in the radiation trans- 
fer process, and the energy transport by conduction and that by radiation 
are independent. For large values of x Q , the radiation transport becomes 
similar to a diffusion process (Rosseland approximation) and* again the inter- 

action with conduction diminishes. 


For the discussion of the effect of gas flow on radiation without con- 
duction it is again desirable to present the results in the form of dimen- 
sionless parameters. In the discussion to follow both plates are at the same, 
constant temperature T* and the pis enters the channel from one end through 
the black porous plug at a temperature Tj_. The gas flow is in the direction 
of length L, as shown in figure 1. The results to be given are for slug 
flow, G is constant. In a manner similar to that described earlier, the 
solution to equation (1) is a function of the dimensionless inlet tempera- 


ture (T^/T # ) and the following set of dimensionless parameters 


The parameter 


£££ 

oT^ 



is known as the Boltzmann number N Bo and was also 


used by Adrianov (3) in presenting his results. In figure 5, the parameter 


T - T. 

A = - 2. , ^ . 1 represents the ratio of the actual heat transferred to the gas 
* ~ i 

to the maximum theoretically possible. The striking feature of figure 4 
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is that for a given value of N-g Q the heat transfer to the gas goes 
through a maximum with increasing t q . The reason for this behavior is 
that as t q increases the layers of gas next to the wall effectively 
absorb most of the radiation from the wall and thereby shield the bulk of 
the gas stream from any direct radiation. Because of this the bulk of the 
gas absorbs less radiation with increasing t q , resulting in reduced 
overall heat transfer. 

The effect of this self shielding can also be seen by comparing the 
transverse temperature profiles at the channel exit for a weak and a 
strongly absorbing gas as shown in figure 6. At t q =1.0 there is very 
little self shielding and the temperature profile across the channel is 
fairly flat. For a higher degree of gas opacity at the same value of A, 
t q = 4. 0, the temperature profile shows the variation expected, the tem- 
peratures next to the wall being higher but the temperatures in the re- 
mainder of the stream being reduced due to the heavy self shielding. The 
preceding remarks have indicated some of the effects of conduction and gas 
flow on radiant heat transfer to absorbing gases. The interaction between 
conduction and radiation resulLs in an increase m net heat transfer in 
comparison to the sum of the conduction and radiation heat transfer taken 
separately. For the case of radiation to a flowing gas from a constant 
temperature surface, the net heat transferred to the gas goes through a 
maximum as the gas absorbtivity is increased. 
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